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ABSTRACT 

I consider the effect of repeating burst sources on the (V/V maa; ) statistic. I 
find that if the burst sources are distributed homogeneously in a <i-dimensional 
space, (V/V m ax) converges to d/(d + 3) as long as the luminosity function 
is independent of position. Choosing the brightest event from a cluster of 
k events shifts the luminosity function to higher luminosities, but if the 
original luminosity function is position- independent, it remains so. Therefore 
the treatment of repeating events, if applied consistently, will not affect the 
effectiveness of (~V/V max ) as a test of burst homogeneity. The calculation of 
(V/V maa; ) for apparent repeating and nonrepeating source populations will be 
biased by the incorrect classification of faint bursts. In conclusion, the current 
practice of calculating (V/V max ) using all bursts, even apparent repeaters, and 
treating multi-spike bursts as single bursts, is valid and consistent. 



Subject headings: gamma rays: bursts 



1. Introduction 



Before the release of the first gamma ray burst catalog (Fishman et al. 1993) from 
the Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray 
Observatory (GRO), the traditional view had been that the "classical" burst sources did 
not repeat while the three known repeating gamma ray burst sources constituted a different 
class, the soft gamma-ray burst repeaters (SGRs), with different spectral properties (see 
e.g., Higdon & Lingenfelter 1990). Recent studies of this catalog suggest that the classical 
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bursts do indeed repeat: Wang & Lingenfelter (1993a,b) have identified a repeating source; 
and Quashnock & Lamb (1993) find an excess of nearest-neighbors at small separation, 
suggesting a population of repeaters. However, Hartmann et al. (1993) do not find the 
expected signatures for repeaters in the angular correlation function or the nearest-neighbor 
statistic. Similarly Narayan & Piran (1993) find only a weak nearest-neighbor excess at 
small separation, and a comparable excess of antipodal bursts (i.e., farthest neighbors 
with separations close to 180°), which they find symptomatic of systematic effects and not 
of repeaters; Maoz (1993) discusses a possible systematic effect. Nowak (1993) identifies 
weaknesses in the nearest-neighbor methodology. Therefore, the existence of repeaters is 
controversial. In addition, Lingenfelter, Wang & Higdon (1993) suggest that complex, 
multi-spike bursts should be treated as a series of repeating events. 

In light of the possibility of repeating burst sources, and the suggestion that each spike 
within a burst should be considered as a burst in its own right, the effect of repeating 
bursts on the (V /V ma x) test for source density uniformity (Schmidt 1968; Schmidt, Higdon 
& Hueter 1988) needs to be evaluated. In this study I use "homogeneous" as uniform in flat 
Euclidean space. I consider a population of events drawn from the same luminosity function 
without regard to the prior history (i.e., each event is independent). Enough sources are 
represented so that space is sampled sufficiently Throughout I deal with the effects of the 
distance to the burst source and not the projected position on the sky. The events from 
each source are clustered, and {V/Vmax) can be calculated either by including all events, or 
only the brightest. 

The cluster of events might be repeated bursts from the same source (i.e., events 
separated by a few hours); since rarely can we identify repetitions, (V/Vmax) is usually 
calculated using each event. The assumption that each burst is independent (i.e., 
characterized by a single luminosity function) is a reasonable assumption in the absence of 
further information. On the other hand, a cluster of events within a few minutes, which are 
undoubtedly from one source given the overall burst frequency, is currently considered to 
be one burst, and (V/V ma x) is calculated using the brightest spike within the cluster. The 
hard-to-soft evolution seen in many bursts (Norris et al. 1986; Band et al. 1992) shows that 
the burst source retains memory of preceding spikes. Therefore the assumption that the 
events are drawn independently from a single luminosity function is mediocre. Nonetheless, 
I assume a single luminosity function for the purpose of the current analysis. 

My argument is as follows: after developing the formalism for calculating (V/V max ) 
(§2), I show that (V/V ma x) is independent of source luminosity functions which are 
functions of luminosity alone (i.e., are not distance-dependent) if the source population 
is uniform in some dimension (e.g., a sphere or a disk). Considering all the events from 
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a repeater, or retaining only the brightest, affects the luminosity function but maintains 
its spatial independence (§3). Separating bursts into apparent repeaters and nonrepeaters 
introduces spatial-dependence into the luminosity function, with (V/V maa; ) biased to lower 
values for apparent repeaters, and to higher values for nonrepeaters (§4). Finally, I draw 
the relevant operational conclusions (§5). 



2. Functional Dependence of V/V r 



Assume bursts have the normalized differential luminosity distribution <f>[L, r] where 
the peak photon luminosity L is in photons-sec -1 , and the burst sources are distributed as 
n[r]. Note that I permit the luminosity function to change with radius. The observed peak 
count rate is C max = L/A-nr 2 . At any given time the detector has a count rate threshold 
Cmin with a normalized distribution g[C min ]. Define £ =C max /C min and v = £ -3 ^ 2 (this is 
the transformation that removes the detector threshold to answer specific questions; see 
Band 1993a). Then the distribution of £, S[£] = dN/d£, is 



m = j ' n[r]g[C min }<f>[L,r}6 
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In brief, the integral in eqn. (1) considers all the sources which can contribute bursts with 
a value £: as the threshold C m i n varies (described by g[C m i n }), different values of C max 
contribute to S; a given C max fixes the ratio of L to r 2 ; the luminosity function <f>[L, r] 
therefore weights the distribution of relevant radii; and finally, the source density n[r] 
provides the number of burst sources at each radius. To get the average V/V maa; =f = £~ 3 ^ 2 
I need the distribution of v 
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where I use dN = Ed!; = \dv (and I drop the minus sign). The average V /V ' max is 

(V/Vma X }= f'dvvxlv}/ f'dvxlv] . (3) 

JO JO 
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Assume n(r) = no(r/r ) a ] clearly the homogeneous 3- dimensional case corresponds to 
a = 0. Then (inserting n(r) into eqn. [2]) 



«/3 f/L\^ ( LV^ \ V2 ' 
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g[C min }dLdC min . (4) 



If the luminosity function 4>[L,r] is a function of luminosity alone, <f>[L], and has no spatial 
dependence then the dependence of x[ v ] on v—V/V max is decoupled from the luminosity 
function (i.e., in eqn. [4] the integrand has no v dependence). In this case 

(v/v maa) = \f; =g±| . (5 ) 

For source distributions homogeneous in a rf-dimensional space n(r)dV oc r d ~ x dr. The 
equivalent 3-dimensional source density which can be used in the above equations for x[v] is 
n(r) oc r d ~ 3 or a = d — 3; this is the density the spatial distribution would have if averaged 
spherically. Consequently, if the luminosity function is spatially independent then 

(V/V max ) = ^3 (6) 

which gives (V/V max )=0.5, 0.4 and 0.25 for d =3 (sphere), 2 (disk — Galactic plane?) and 1 
(linear — spiral arm?). An important corollary is that any manipulation that maintains the 
spatial independence of the luminosity function will not affect (V/V mox ). 

In addition to its average values for homogeneous source densities, V/V max will be 
distributed uniformly between and 1 for d — 3, providing a secondary test for this 
dimension. Variants tailored to dimensions other than d — 3, such as the (A/A max ) test for 
d = 2 (Kluzniak 1992), can be constructed by transforming to £ _d//2 . 



3. The Luminosity Function of Burst Clusters 



Assume repeating sources all have the same spatially independent luminosity function 
4>i(L), and that the source sample is large enough so that true averages are observed. The 
repeating sources can be treated either by considering only the brightest event in calculating 
(V/V mox ), or by including all observed events. I will show that the luminosity function for 
each case is spatially independent, and will find the effect these choices have on (V/V max ). 

First I find the luminosity function 4>k[L] for the brightest event in a cluster of k 
events. Since the differential luminosity function is normalized, the cumulative luminosity 
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function $i[L] varies uniformly between and 1 (by definition). Therefore this cumulative 
distribution is a mapping of the luminosity into a uniform distribution: 

gAu] = 1 for < u < 1 where u= fa[L'\dL' = fcJLl . (7) 

Jo 

If there are fc events in a cluster, then the distribution of the value of the single event 
cumulative luminosity function corresponding to the brightest event is 

g k [u] = ku k - x for < u < 1 (8) 

because there are k events which could be the brightest and there are k — 1 other events 
which must be less bright; the probability that an event is less bright is u. Note that since 
the average value of u is 

(«) = J 9k(u)udu = , (9) 

the average luminosity is shifted to higher values, as expected. The increase is not very 
great for flat luminosity functions (e.g., by no more than a factor of 4 for fa oc L^ 1 / 2 ), but 
can be very dramatic for steep luminosity functions (e.g., by more than a factor of 50 for a 
ten-event cluster and fa oc L~ 2 ). The luminosity function of the brightest burst is derived 
by relating it to the single event luminosity function. The cumulative luminosity function 
of the brightest event is also distributed uniformly: 

G[U] = 1 for 0<U<1 where U = [ L fa\L']dL' = $JL] . (10) 

J o 

Equating G[U]dU = gk[u]du gives 

U = u k 

J q fa\L']dL' = U fa[L']dL'\ (11) 
fa[L] = fc$i[L] fc "Vi[£] . 

Note that is the luminosity function regardless of whether all k events can be detected. 
Similarly, since each event is assumed to be independent, a cluster of M events can be 
considered as two clusters of N and M — N events if the separation into two clusters is 
based on a luminosity-independent criterion, e.g., the temporal separation between events. 

Let fk be the fraction of clusters with k events; obviously Y^k=i /fc = 1, and the average 
number of events J2T=i k fk must converge. If each event is considered separately, 

oo 

<f>[L]=fa[L}J2kf k (12) 
k=i 
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where the average number of events has been inserted to preserve the total number of 
events observed (n[r] is the density of sources). However, if we consider only the brightest 
event in each cluster, then 

oo 

4>[L] = £ fML] ■ (13) 

k=l 

Eqn. (13) is applicable if the events from a repeating source are broken into a number of 
clusters based on a luminosity-independent criteria. For example, events separated by less 
than some predetermined time (e.g., half an hour) may be considered a single burst, events 
after a gap greater than this time are considered as part of a new burst. 



4. Separating Apparent Repeaters and Nonrepeaters 



There is a temptation to separate repeating burst sources from nonrepeaters, and 
calculate (V/V max ) for each class, perhaps to determine whether they originate from 
different source populations. However, repeaters can only be identified if the second 
brightest event exceeds the detector threshold. Whether the repeating source is identified as 
such therefore depends on the distance to the source. This introduces a spatial dependence 
into the luminosity functions, complicating the determination of the source distribution. 
Sources observed to repeat will be closer, and have a smaller (V/V mox ). 

Consider a cluster of k events with the luminosity of the brightest event corresponding 
to a value u\ of the single event cumulative luminosity function (see eqn. [7]). Then u 2 , 
the value of the single event cumulative distribution for the second brightest burst, is 
distributed as 

9kA u 2 I «i] — (k — l)i4~ 2 u\~ k , u 2 <U! ; (14) 

there are k — 1 possible second brightest events, and the k — 2 less bright events must have 
values of u less than u 2 . The dependence on u\ results from restricting all k — 1 bursts to 
the range to ui. Then the fraction of the clusters for which the brightest event has a value 
Mi and the second event has a luminosity which produces a count rate less than threshold 
(i.e., L 2 < L = Anr 2 C min ) is 

f s ,k = / (k-l)vt 2 u\- k du 2 = (u /u 1 ) k - 1 
Jo 

\l L ML'}dL'J { JoML'}dL> ) 
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where u = $i[L ] and L = 4irr 2 C min = ^ir^Cmax^ 1 = Lv 2 ^. Thus / Sjfc is the fraction 
of the /c-event clusters which are counted as nonrepeaters (i.e., misidentified) because the 
second spike is below threshold. Note that f s ^ is a function of both the intensity of the 
brightest event, and the distance to the source (through Lo). The effective differential 
luminosity functions for apparent single events <p s and for apparent multiple events (f> m are 

oo oo 

4>s = h<t>i + fs,kfk<t>k and <f) m = ^2(l- f Sik )fk(f>k ■ (16) 

k=2 k=2 

Because f s ^ is a function of distance to the burst source, the luminosity function for the 
apparent repeaters and nonrepeaters is also spatially-dependent. Since f s ^ increases with v, 
it causes k-event clusters to be misclassified as apparently single events with larger values 
of (V/V max ) and correctly classified as repeaters with smaller (V/V max ). 



5. Conclusions 

The {V/V m ax) statistic is best used to test whether sources are distributed uniformly in 
Euclidean, 3-dimensional space (Band 1993a), or by extension, in any other <i-dimensional 
space. I showed that this diagnostic power assumes the luminosity function is spatially 
independent. Therefore, transformations which do not introduce a spatial dependence into 
the luminosity function will not affect the utility of (V/V max ). 

Using only the brightest event from a cluster of k events modifies the effective 
luminosity function but maintains its spatial independence. As expected, the average 
luminosity increases with the number of events, which means that (V/Vmax) probes greater 
distances, and is therefore more sensitive to inhomogeneities in the source distribution. 
However, if the source distribution is indeed homogeneous, the treatment of event clusters 
does not affect {V/V max ). 

Apparent repeater and nonrepeater populations do not have meaningful values of 
(V/Vmax) ; (V/Vmax) is biased towards smaller values for apparent repeaters and towards 
larger values for apparent nonrepeaters. The second brightest event from a repeater must be 
detectable, and therefore the fraction of true repeaters which are classified as nonrepeaters 
increases with distance. Thus multi-spiked bursts should not be separated from single-spike 
bursts, and repeating sources from nonrepeating. 

The burst phenomenon may consist of events with a distribution of separation times 
peaked at short timescales, but with a tail to large separations (Lingenfelter et al. 1993). 
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Those events within a few minutes of each other are treated as a single burst, while those 
separated by hours or more are considered different bursts. Note that the existence of 
repeating bursts separated by more than a few minutes is controversial, as discussed in 
the Introduction. This separation into individual multi-spike bursts and repeating bursts 
does not depend on the distance to the burst source, and should not introduce a spatial 
dependence into the luminosity function. Therefore the current method of calculating 
(V/V ma:r ) does not compromise the test for homogeneity. 

There are two alternatives to the current methodology: calculating (V/Vmax) using 
each spike in a burst or including only the strongest burst from a repeating source. If 
performed consistently, each alternative will not affect the diagnostic power of (V/V mox ). 
However, the assumption that in a multi-spike burst each spike is drawn independently from 
the same luminosity function is invalidated by the observed spectral evolution. Similarly, it 
is probably impossible to identify distant repeating sources or faint repetitions from nearby 
sources as a consequence of the large uncertainty in the projected position on the sky of low 
count rate bursts; incomplete removal of all but the strongest event from a repeater will 
bias (V/V m ax) to higher values. Therefore these variants do not increase the effectiveness 
of (V/Vmax), and can easily introduce biases. 

Previously I (and others) showed that (V/Vmax) is appropriate for testing the 
hypothesis of source homogeneity, but not for investigating inhomogeneous source 
distributions (Hartmann & The 1992; Petrosian 1993; Band 1993a). Maximum-likelihood 
methods using the observed peak flux and the known detection threshold distributions 
(Band 1993b; Loredo & Wasserman 1993), survival analysis (Efron & Petrosian 1993), 
and moment methods (Horack et al. 1993) are more powerful for studying inhomogeneous 
sources. Since BATSE has found that (V/V ma x) is definitely not 1/2, and the intensity 
distribution is not a power law (e.g., Meegan et al. 1993; Fishman et al. 1993), the burst 
source distribution is clearly not homogeneous. Therefore, (V/V ma x) should be retired 
in favor of methods which utilize all the information in the burst distribution. For these 
superior methods the treatment of repeaters discussed here is just as relevant. Since the 
luminosity function is unknown and must be modeled, manipulations which alter the 
luminosity function will not complicate significantly an already difficult problem. However, 
it is preferable to keep the luminosity function distance- independent. 

In summary, the current practice of treating all events which are separated by less 
than a reasonable time (e.g., ~half-an-hour) as one burst and retaining all the bursts from 
an apparent repeater produces a burst database appropriate both for the calculation of 
(V/Vmax) and for studies of inhomogeneous source models. 
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